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Tamami Kusuda and Katsumi Ishii 


ABSTRACT 


The majority of the available solar radiation data for the United States and Canada are 
in the form of either monthly averages of daily total radiation on a horizontal surface or 
hourly values for cloudless days on horizontal, tilted, and vertical surfaces. Hourly solar 
radiation data for walls and roofs under "average" solar conditions were computed in order to 
be able to make estimates of the effect that incident solar radiation has on a building and/or 
its heating and air conditioning system over a heating and/or cooling season. Calculation 
procedures developed by Liu and Jordan were modified and hourly solar radiation data were com- 
piled for 80 locations in the United States and Canada. Data were also computed and tabulated 
for a new parameter called "sol-air temperature for glass". 


Key Words: Computer; energy; radiation; solar. 
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Surfaces on Average Days in the United States and Canada 


by 


Tamami Kusuda and Katsumi Ishii 


1. Introduction 


There are basically two ways that solar radiation data can be used for design purposes 
in the building profession. One way is to use past hourly data for the locality or similar 
locality in question, estimate through mathematical simulation what the performance factors 
of the building and for its components would have been under those conditions, and base future 
predictions for performance on that analysis. The second way is to use data that have been 
derived for specific types of days, either completely clear, "average", or statistically ana- 
lyzed data for a particular solar process, to calculate the performance of the building and 
for its components for those types of days. 


The first technique has the advantage of easily accounting for the detailed transient 
nature of the interaction between the solar radiation and the building, and the building and 
its heating and air conditioning system. This is particularly important when the heating and 
air conditioning system has been specifically designed to collect the solar energy and utilize 
it inside the building for heating and/or cooling. ASHRAE has been very active in the recent 
past in outlining recommended procedures to conduct such hourly calculations for computing 
heating and cooling loads on buildings [1]*. In the ASHRAE-recommended procedures [1], com- 
pletely clear day data are predicted using theoretical equations and these data are in turn 
modified with actual hourly recorded cloud-cover data. Computer programs that have been: 
written for predicting the performance of solar heating and/or cooling systems on an hourly 
basis use actual recorded hourly solar radiation data [2, 3]. 


The primary disadvantages of using this first technique are the cost and complexity asso- 
ciated with the use of the computer programs and the lack of good hourly data, particularly 
for the solar radiation itself. 


The major source of solar radiation data in the United States is the National Weather 
Service and its network of stations. The data collection began in 1902 with the first mea- 
surements of the direct component of the solar radiation at normal incidence being made in 
Asheville, North Carolina. At the present’ time, there are 86 stations covering the U.S. and 
the West Indies. Reference 4 gives a complete history of the network. The majority of the 
present 86 stations record only the total radiation on a horizontal surface. At five of the 
stations, the direct component is also measured. Of the 86 stations, 6/7 have their data pro- 
cessed for daily total only, while hourly data are processed at the other 19. Reference 5 
describes how the data can be made available either on tape or punched cards. 


In addition to the relatively small amount of hourly solar radiation data generally avail- 
able as described above, the quality of the data has been questioned. Jessup has reported [4] 
that estimated data errors range from + 5 percent to + 30 percent. This has been mainly due 
to insufficient resources in ‘all but the recent past to properly maintain the network. A com- 
plete analysis of some of the specific problems with the network instrumentation is included 
in references 6 and 7. 


There is evidence, with the increasing emphasis on the use of solar energy in buildings, 
that much of the past solar radiation data will be rehabilitated and additional data will be 
collected in the future; however, it is unlikely that the hourly data to be taken will be ex- 
tended to cover surfaces other than the horizontal for the majority of the stations. For this 
and other reasons given above, emphasis will be placed in the foreseeable future on using solar 
radiation data for specific types of days and orientations for design purposes. 


See references at end of text. 


The 1977 ASHRAE Handbook of Fundamentals [8] will include tables of solar heat gain fac- 
tors for all combinations of six different northern latitudes, sixteen primary compass orien- 
tations plus the horizontal, for the 2lst day of each month of the year and for each hour of 
that day. This enables one to calculate the amount of solar energy passing through a unit 
area of reference glazing (single sheet of 0.3 cm (1/8 in.) thick double-strength window glass) 
on a completely clear day, and is based on a calculation procedure developed by Stephenson 
[9]. Recently Morrison and Farber [10] have extended this type of calculation and have pro- 
duced tables of incident solar radiation for all combinations of six different northern lati- 
tudes, five different south-facing surfaces tilted up at angles ranging from latitude minus 
10 degrees to the vertical plus the horizontal, for the 2lst day of each month of the year 
and for each hour of that day. Both of these above sets of data can be used for design pur- 
poses to determine the maximum effect that the solar radiation is likely to have for a given 
condition or in a given situation. However, it is felt desirable to also have similar tables 
of data for "average" conditions so that estimates can be made for the long-term effect that 
incident solar radiation has on a building and for its heating and air conditioning system. 


References 11 and 12 contain a series of 12 monthly maps showing isoclines of daily total 
radiation on a horizontal surface for the United States and the world, respectively. The data 
were gathered from a large number of sources and can be used to indicate general trends and 
to predict areas of potential solar energy applications. 


Liu and Jordan [13, 14, 15] conducted extensive analyses during the early 1960's on the 
available solar radiation data and developed several empirical correlations that can be used 
to estimate the available solar radiation on "average" days for each month of the year and 
for a large number of locations in the United States and Canada. Using the correlations, it 
is possible to take the monthly average daily total radiation on a horizontal surface, divide 
the daily total into direct and diffuse components, convert each component into hourly values, 
and then compute the hourly value of either component on a surface of any orientation desired. 


This report outlines the technique that was used to compute and tabulate the monthly 
"average" incident radiation on an hourly basis during the day for each month of the year, 
and each of eight different vertical orientations. The data was tabulated for 80 different 
locations in the United States and Canada. A minor modification was made to the basic tech- 
nique recommended in the Liu, Jordan procedure [14, 15] in order to eliminate what was believed 
to be erroneous results near the sunrise and sunset hours. An additional parameter called 
"sol-air temperature for glass'' was also computed and tabulated for each of the locations. 
An explanation of the meaning of this parameter is included in Section 7 of the report. 


2. Calculation Procedure of Liu and Jordan 
The pertinent correlations developed by Liu and Jordan [13, 14, 15] were for: 


a. the monthly average daily diffuse radiation on a horizontal surface as a func- 
tion of the monthly average daily total radiation on a horizontal surface, 


b, the ratio of the hourly diffuse radiation to the daily diffuse radiation (both 
on a horizontal surface) as a function of the number of daylight hours, and 


c. the ratio of the hourly total radiation to the daily total radiation (both on 
a horizontal surface) as a function of the number of daylight hours. 


The correlations were based on available solar radiation data from Blue Hill, Massachusetts; 
Nice, France; Helsingfors, Finland; and Kew Observatory, London, England. 


Liu and Jordan [15, 16] also compiled for each of twelve months at 80 different locations 
in the United States and Canada, the "average" daily total radiation on a horizontal surface 
(H) and the ratio Ky = H/Hpo, where Hp is the daily total extraterrestrial radiation on a hori- 
zontal surface calculated for the 15th or 16th of each month. The ratio Kr represents the 
fraction of extraterrestrial radiation that is transmitted through the atmosphere and thus 
provides a useful index for the solar climate of a particular location. 


The correlations described in a, b, and c above are included herein as Table Lig Figure 
1, and Figure 2, respectively. Ka in Table 2 is defined as D/H, where D is the "average"! 
degigy diffuse radiation on a horizontal surface. The step-by- Been technique for using the 
correlations to compute the "average" incident solar radiation on a horizontal or vertical 
surface is as follows: 


a. Daily total extraterrestrial radiation on a horizontal surface for a given 
latitude and declination angle is calculated as follows: 


H =—riI (cos L cos §S sin W +W_ sin L sin S) (1) 
T sc s 8 


where 


r = ratio of solar radiation intensity at normal incidence outside 
the earth's atmosphere to the solar constant (monthly average 
values are given in Table 1) 


I__ = solar constant, 1,393 W/m? (442 Btu/(h + ft?) * 
L = latitude, radians 


S = solar declination angle, radians, monthly values are given in 
Table 1 


W_ = sunrise hour angle, radians, obtained from: 


cos We = - tan L tan S (2) 


b. Obtain the monthly average daily total radiation on a horizontal surface, H, 
from references 15 and 16 or other available sources. 


ce. Calculate K by: 
a oe 
KF (3) 
) 


or obtain the value from references 15 and 16 for any of 80 locations. 


d. Determine the value of K, from Table 2. 


e. Calculate the monthly average daily diffuse radiation on a horizontal surface, 
D, by: 
D= 6 4 
DeaRken AHy (4) 


f£. Calculate the value of a factor, Ty from the following equation or determine 
from Figure 1: 


cos W - cos We 


d ¥ 24 sin W.-W. cos W (5) 
s s s 


r 


where 


= hour angle corresponding to a given hour of the day (solar time), 
radians 


* 
This was the value used by Liu and Jordan. A value of 1,353 W/m (428 Btu/(h - £t7)) is 


the value currently recommended by a number of governmental agencies and technical soci- 
eties in the U.S. [17]. 


Calculate the average hourly diffuse radiation on a horizontal surface a) 
by: 


Ty Sir ened); (6) 


Determine the value of a factor, tp from Figure 2. 


Calculate the average hourly total radiation on a horizontal surface ( Th) 


by: 


I =r H (7) 


Calculate the average hourly direct radiation on a horizontal surface Ge 


) 
by: n 


=) ee (8) 


Calculate the angle of incidence between the sun's direct beam and a hori- 
zontal (84) or vertical surface (8) by: 


cos OF = cos L cos S cos W + sin L sin S (9) 
or 
cos oF = cos ALT cos SAZM (10) 
where 


ALT = solar altitude angle, the angle between the sun's direct beam 
and the horizontal, radians 


SAZM the difference between the solar azimuth angle (AZM) and the 


wall azimuth angle (WAZ), radians 
Figure 3 shows the relationship between the various angles. 


Calculate the average hourly direct radiation on a given vertical surface 


(1 by: 


I 
= Dh 
Ly = cos 6, + Cos oF (11) 


Approximate the average hourly diffuse radiation on a given vertical surface 
(E>); by: 
dv 


ap 
iy aude 
hee (12) 


Approximate the average hourly reflected radiation from the ground on a given 
vertical surface Coe) by: 


Tre _£ Th (13) 
where 


ie = ground reflectance 


o. Calculate the average hourly total radiation incident on a given vertical 
surface by: 


mT seit. a (14) 


3. Anomalous Calculated Values Using the Liu and Jordan Calculation Procedure 


The procedure described in the previous section was used to generate the hourly solar 
radiation values incident on vertical surfaces oriented in the eight primary directions (S, 
SW, W, NW, N, NE, E, and SE) as well as the horizontal surface. Upon analyzing the data, 
exceptionally large values of incident solar radiation were observed near the sunrise and 
sunset hours. A close examination of the calculation procedure revealed that the large peak 
resulted from the use of equation (11) where cos 6y is very close to zero at the sunrise and 
sunset hours. This phenomenon is illustrated in Figure 3 where the value of the direct com- 
ponent of the solar radiation (Ipy) is computed for a vertical surface oriented directly fac- 
ing the sun (equation (11) with cos e., = 1) on February 22 in Miami, Florida: 


i In = Tyn/ 698 On (15) 


i The inconsistency appeared to arise as a result of the sensitive empirical relationship between 
Kr, rp, Kg, and rg. Upon discussions with one of the authors of [15, 16]*, it is felt that 

the inconsistency is mainly due to the inaccuracy of the empirical relationship for rp (Fig- 
ure 2). 


4. ASHRAE Methodology for the Computation of the 
Direct Component of the Solar Radiation 


In order to circumvent the anomaly described in the previous section, it was decided to 
replace equation (15) with an equation similar to one used in reference 18: 


utr Se) 
| ee ath cos On 

= 16 
Coe (16) 


The original equation was developed by Stephenson [9] for completely clear sky conditions. 
It has been assumed in this work that the same basic form of the equation is valid for 
"average" day solar radiation. A modified multiplying coefficient has been adapted, denoted 


in equation (16) as A (as opposed to A in reference 9). 


The monthly values of B for equation (16) are given in Table 1. In order to develop 
appropriate values of A, it_was necessary to calculate the average daily direct radiation 
over a horizontal surface, Hy by: 


sunset 
AD, = Hoe cos On dt (17) 
sunrise 


where t denotes time. Combining equations (16) and (17): 


x 
Private communication with Dr. Liu. 


W 
s 
al OV ie 
Hoh * ae A J e cos en dw (18) 
0 


By denoting the integral of the above equations by G, and by using K and K, of the Liu and 
Jordan procedure, equation (18) can be written: 


el aie 
Hh A G2 & - KE) (19) 


By substituting equation (1) into equation (19), the value of A can be calculated as follows: 


A= (Ky - Ky) f (20) 
where 
r I. 
f= . (cos L cos § sin Ws + We sin L sin S) (21) 


Table 3 lists the values of f which were computed using the above equations. The values 
are tabulated for latitudes of 0 through 50 degrees in one-degree increments and for each of 
the 12 months. The monthly values are plotted in Figure 5. 


The concept of f is convenient since its value is independent of the specific local cli- 
mate. As an example, both Hatteras, North Carolina and Albuquerque, New Mexico are at a lat- 
itude of 35° N. Table 3 indicates the value of f in January to be 644.5 at this latitude. 
The Liu and Jordan tabulated data [15, 16] show the following values of Kh and Ka for these 
two cities. 


Hatteras Albuquerque 
K 0.634 0.704 
K, 0.166 0.147 


Using these data, A for Hatteras and Albuquerque were determined as 302 and 359 respectively, 
whereas the cloudless sky value of A developed by Stephenson [9] was 390. 


5. Modification of the Liu and Jordan Procedure to Eliminate 
the Anomalous Calculated Solar Radiation Values 


be With the relationship between equation (16), the constant A, and the constants Ke and 
Ky established, a modifted hourly calculation procedure was developed as follows: 


a. Same as before. 
b. Same as before. 
c. Same as before. 
d. Same as before. 


e. Same as before. 


f. Same as before. 

g- Same as before. 

h. Determine the value of f from Table 3 for a given month and latitude. 
i. Calculate A using equation (20). 


j. Calculate the hourly value of ION using equation (16). 
k. Same as before. 


1. Calculate the average hourly direct radiation on a horizontal surface by: 


Toh = Th cos 04 (22) 


m. Calculate the average hourly total radiation on a horizontal surface by: 


Sr es (23) 


n. Calculate the average hourly direct radiation on a vertical surface by: 


Ha Ion cos 6. (24) 


o. Same as m. in Section 2. 
p. Same as n. in Section 2. 


q- Same as o. in Section 2. 


6. Computation and Tabulation of Data for Eighty Locations 


A computer program was written for carrying out the calculation procedure described in 
the previous section. A listing of the program is given at the end of this report. Complete 
tables of total solar radiation values are also included at the end of the report in Table 4 
for 80 different locations in the United States and Canada [15, 16]. The tabulated solar radi- 
ation data are all in units of either Btu/(h + ft%) (= 3.152 W/m2) or Btu/ft? (= 11,348.9 J/m2). 


Sample results are shown in Figures 6 through 12 for the horizontal and different verti- 
cal orientations in February in Miami, Florida. Shown in these figures are the hourly radi- 
ation values calculated by Liu and Jordan, the authors' modification of the Liu and Jordan 
procedure, and the data for completely clear skies [10]. 


7. Calculation of "Sol-Air Temperature for Glass" 
In general, energy is transferred through glazing both by transmitted solar radiation and 
by conduction and convection as a result of a temperature difference between the outside and 


the inside air. The governing equation for the net heat transfer, q (W), across the glazing 
into a building space is: 


CONG RS UIG Sey] (25) 


where 
A = cross-sectional area of the glazing, m” 


effective transmittance of the glazing to short-wave radiation 


a| 
ul 


I = incident total solar radiation on the glazing, W/m? 

U = overall heat transfer coefficient for the glazing, (W/ (m2 * °C)) 
t_ = outside air temperature, °C 

t, = inside air temperature, “AG 


If the first term in equation (25) is multiplied and divided by U, the equation can be re- 
arranged to: 


Tal 
Q=UA | 7 ate team “| (26) 
or 
q=UA [sare = S| (27) 
where 
et tae 
SATG = 7 = sol-air temperature for glass, °C 


Values of SATG are useful in assessing the net heat transfer through glazing into the build- 


ing. In other words, for those periods of time when SATG > t;, there is a net heat gain into 
the building and when SATG < t,, there is a net heat loss. 


Equation (26) is written on an instantaneous basis and can be used for computing hourly 
values. The equation can also be written to be valid for a period covering a number of hours 
N: 


aN, 
Teds a 
Q =aueal| ===" Na eaNene) (28) 
U i 
where 
Q = net heat transfer over the time period N, W 
N = number of hours 
ae = average outdoor air temperature over the time period N, °C 


The indoor air temperature, ty, has been assumed to be constant for the t‘me period N. Equa- 
tion (28) can be rearranged as was done with equation (25): : 


aN 
Taras IE 

Q=UAN +7 -t, (29) 

Q=UAN [sac - =| (30) 


where 
is 
qe I 
i 


SAG = ape A&C (31) 


ee 2h a 
UN fo) 


The values of SAG were computed for each month of the year, for each of the eight primary ver- 
tical orientations and the horizontal, and for each of the 80 locations in the United States 
and Canada. The values are tabulated in Table 4 in units of °F along with the average hourly 
radiation values discussed previously. 


Four different values of SAG were computed. Two values were for a 24-hour period (N = 
24), one for single glazing (SATGS) and a second for double glazing (SATGD). The other two 
computations were also for single and double glazing, but the time ie eee N was shortened to 


include only the daylight hours. In all cases, the radiation values rel I were the total inci- 


dent radiation under "average" conditions discussed previously. The values of t used in the 
computations were 0.87 for single glazing and 0.78 for double glazing. The overall heat trans- 
fer coefficients were 6.4 W/(m2 + °C) (1.13 Btu/(h - ft2 + °F)) and 6.0 W/(m2 + °C) (1.06 Btu/ 
(h . ft2 - °F) for single glazing and for winter and summer, respectively. The values used 

for double glazing were 3.1 W/(m2 - °C) (0.55 Btu/(h + ft2 + °F)) and 3.1 W/(m2 + °C) (0.54 
Btu/(h - ft2 - °F), respectively. The value of tg was either the average over the 24-hour 
period or daylight hours only, depending on the time period N. 


The intent in calculating and compiling the data was to permit computation, for a given 
glazing size and orientation, of the net heat transfer into the space of a building for 
"average" solar conditions and for any period of months desired. By repeating the calculation 
for all the glazing in a given building, the total net effect for the entire building could be 
determined. 


8. Summary 


The Liu and Jordan method for calculating the incident hourly solar radiation on a sur- 
face of any tilt and orientation under "average" conditions was modified and used to compute 
the average hourly incident total radiation on the eight primary vertical orientations and 
the horizontal for each month of year and for each of 80 locations in the United States and 
Canada. In addition, a new parameter was introduced and computed for the 80 locations called 
"sol-air temperature for glass''. The tabulated values enable one to easily compute the net 
heat transfer through single or double glazing under "average" solar conditions for any of the 
eight primary vertical orientations or on a horizontal surface. 
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Figure 1 Ratio of Hourly Diffuse Radiation to Daily Diffuse Radiation 
as a Function of Time of Day [15, 16] 
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Figure 2 Ratio of Hourly Total Radiation to Daily Total 
Radiation as a Function of Time of Day [15, 16] 
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Figure 4 Anomalous Sunrise Hour Solar Radiation Generated 
by the Liu and Jordan Procedure for February 22 
in Miami, Florida 
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Figure 5 Values of Monthly Average Daily Solar Factor f for Various Latitudes 
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Figure 8 Comparative Plots of Hourly 
Solar Radiation Incident 
Upon South-facing Walls for 
February 22 in Miami, ‘Florida 
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Figure 9 Comparative Plots of Hourly 
Solar Radiation Incident 
Upon North-facing Walls for 
February 22 in Miami, Florida 
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Figure 10 Comparative Plots of Hourly Solar Radiation Incident Upon East- 
(or West)-facing Walls for February 22 in Miami, Florida 
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Figure 11 Comparative Plots of Hourly Solar Radiation Incident Upon Southeast- 
or Southwest-facing Walls for February 22 in Miami, Florida 
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Figure 12 Comparative Plots of Hourly Solar Radiation Incident Upon Northwest- or 
Northeast-facing Walls for February 22 in Miami, Florida 
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Table 1 


Solar Declination, S, the Ratio of Solar Radiation Intensity 


at Normal Incidence Outside Earth's Atmosphere to the Solar 
Constant (vy), and Atmospheric Extinction Coefficient B [18] 


Month 


January 
February 
March 
April 
May 

June 
July 
August 
September 
October 
November 


December 


Declination 


21 


Ratio 


1.0300 


1.0207 


1.0057 


0.9875 


0.9727 


0.9670 


0.9692 


0.9785 


0.9945 


1.0133 


1.0267 


1.0327 


Extinction 
Coefficient 


A 


Table 2 The Relationship Between the Monthly Average Daily Diffuse 
and the Monthly Average Daily Total Radiation [15, 16] 


O.3 0.4 0.5 0.6 0.7 


0.179 0.183 0.188 0.174 0.149 
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Table 3 The Value of Factor (A/ (Ky ~ K,) 


559. 
559. 


560. 
560. 
561. 
561-6 
5626 
563. 
5636 


564. 
5656 
566.6 
567-6 
567-6 
568-6 
5706 
S71. 
5726 
573-6 


574. 
5766 
577.6 
579. 
580. 
582-6 
584.6 
586.6 
588. 
590. 


592-6 
594. 
597. 
599. 
602. 
605-6 
608-6 
611. 
614. 
618. 


622-6 
6266 
630. 
635-6 
640. 
645.6 
651. 
657.6 
664. 
671. 


678. 


3 


558.6 
558-6 
5586 
5586 
558.6 
558.6 
558. 
559. 
559-6 
559. 


559. 
5606 
5606 
5616 
561. 
5626 
5626 
5636 
564. 
564. 


565-6 
5666 
567. 
568.6 
5696 
5706 
S71. 
5726 
5736 
575. 


576-6 
578. 
5796 
581. 
5826 
584. 
5866 
588.6 
590-6 
5926 


595-6 
597. 
600. 
602.6 
605-6 
608. 
6126 
615-6 
6196 
622-6 


6266 


MONTH 
6 


588. 
587.6 
5866 
5856 
584. 
5836 
5826 
581. 
580.6 
5796 


5796 
578. 
578.6 
577. 
577.6 
5766 
5766 
5766 
5756 
5756 


5756 
5756 
5756 
5756 
575.6 
5756 
5756 
5766 
5766 
5766 


S06 
577. 
5786 
578. 
5796 
579-6 
5806 
5B1le 
5826 
5836 


584. 
585-6 
5866 
587.6 
589. 
590-6 
591.6 
5936 
595. 
596.6 


598-6 
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587.6 
5866 
585.6 
584. 
584. 
583. 
58?. 
581. 
5816 
580. 


580. 
579. 
5796 
578.6 
578.6 
578. 
578.6 
578. 
578.6 
577.6 


577.6 
57Re 
578. 
57B.6 
578.6 
578. 
579. 
579. 
579. 
580. 


5806 
581. 
582. 
582. 
583-6 
584. 
585-6 
586-6 
587.6 
588. 


589. 
591. 
5926 
593-6 
5956 
596.6 
598. 
6006 
602. 
604. 


6066 
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5B. 
5806 
5B06 
579.6 
579-6 
5796 
57Be 
57B. 
5786 
578.6 


5786 
5786 
578. 
578.6 
578-6 
5786 
578. 
578.6 
5796 
579.6 


579. 
580-6 
5806 
581. 
581.6 
582. 
583-6 
584. 
5846 
5856 


586.6 
587.6 
588. 
5906 
5916 
5926 
594. 
595-6 
597.6 
598. 


6006 
6026 
604-6 
606-6 
608-6 
610-6 
613-6 
6156 
618. 
6216 


623-6 


5i71e 
5726 


574.6 
5756 
5766 
577.6 
578. 
579. 
580. 


582. 
5836 
584. 
586. 
588. 
589. 
591. 
593. 
595-6 
597. 


599. 
602-6 
604. 
607. 
609. 
612. 
615. 
619. 
622. 
6266 


6306 
634. 
638.6 
643. 
648. 
653-6 
659. 
6656 
6726 
679. 


687-6 
695. 
704. 
714. 
725.6 
737. 
7506 
764. 
780. 
798. 


819. 


5726 
5736 
574 es 
575 s 
576.6 
577. 
578. 
579. 
581. 
582.6 


584. 
585-6 
587. 
588. 
590-6 
592.6 
594. 
5966 
599. 
601. 


604. 
6066 
609. 
612.6 
615. 
619.6 
622.6 
626-6 
6306 
634. 


639-6 
644. 
649. 
655.6 
661. 
667.6 
674. 
682. 
690- 
699. 


708. 
719. 
731.6 
744, 
7586 
774. 
7926 
812. 
835. 
B61. 


891. 


Program for Generating Table 3 


DIMENSION DDEC(12) 7-19-51 9-30.28 9002011 e569 20 01 4923.279 202026912003 


*0900637 9-10 0479-19258 9-23.62 7/oDLAT(S1) oFAC(12) 


DIMENSION B(12)7.1420 61440 01560 01800 61967 02059 02079 2010 01779 1609 


*o1499.142/ 


DIMENSION R(12)/71-0391-020791 200570 698750 097278 09670 » 96920 097859 


* 9945012013391 6026791-2-0327/ 
PI=3.1416 

SC=442.1 

WRITE (6°60) 


60 FORMAT(/10X010("**) +" THE VALUE OF FACTOR ( A/(KT=KD) 


*) 
WRITE (60 70) 
70 FORMAT (/5Xe *LATs*27Xe *MONTH® ) 
WRITE (6080) 
80 FORMAT (/5Xe "DEG 1 2 3 uy 
* 11 12") 
WRITE (6081) 
81 FORMAT () 
DO 10 J=1°51 
DLAT(I)=I-1. 
XLAT=INT(DLAT(I)) 
RLAT=(XLAT+(DLAT( I) —=XLAT)/0-6) *PIT/1A06 
DO 20 J=1912 
XDEC=INT(DDEC (J) ) 
RDEC=(XDEC+ (DDEC (J) =XDEC) /0.6) *PI/130.6 
COSWS==TAN(RLAT) * TAN(RDEC) 
wS=ACOS(COSWS) 
OWS=WS¥*180./PI 
COSLD=COS(RLAT) *COS(RDEC) 
SINLD=SIN(RLAT) *SIN(RDEC) 
S=0.6 
DO 30 N=1*39 
W=WS*N/40 
COSZ=COSLD*COS (WwW) +SINLN 
PAR==B(J)/CO0SZ 
APA=ABS (PAR) 
IF(APA.GTe80) GO TO 4O 
ANS=EXP (PAR) *COSZ 
60 TO 50 
40 ANS=0. 
50 S=ANS+S 
30 CONTINUE 
ANSO=EXP (=B (JU) /(COSLD+SINLD) )*(COSLD+SINLD) 7/2. 
AINT=WS/40-* (ANSO+S) 
THOE=R (J) *SC* (COSLD*SIN (WS) +WS*SINLD) 
FAC(J)=THO/AINT 
20 CONTINUE 
WRITE(6090) OLAT(I) o FAC 
90 FORMAT (5X0eF3.002X%012F5-0) 
DO 1 K=10»50010 
IF(I-EQeK) WRITE (6081) 
4 CONTINUE 
10 CONTINUE 
STOP 
END 
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Program for the Average Day Hourly Solar Tables 


P MAIN PRIGRA™M 


DIMENSION H(12) 62KT(12)2TO(12) eTOT(12) 
DIMENSTON TOWN(4) 


DO 1 K=181 

READ(5+4) TOWN 

READ(5+2) He ZKTrTO*XLATeCITYe TOT 

IG=1 

SC=0.9 

RHO=0.2 

CALL LIU(Z2KT+He TOrIGeSCeXLAT# RHO» TOWN? TOT) 
CONTINUE 

FORMAT (12F6.1/12F6.3/12F60190F 5.2 0A3/12F 401) 
FORMAT (1H1) 

FORMAT (4A6) 

STOP 

END 


SUBROUTINE LIU(ZKT+H» TO» TGeSC*XLAT* RHO» TOWN? TOT) 
THIS SUBROUTINE CALCULATFS MONTHLY AVERAGE SOLAR HEAT GAIN THROUGH WINDOW2 
OF FOUR DIFFERENT ORIENTATIONS. 
ZKTeeeLIU/JORDAN FACTOR = DAILY TOTAL RADIATION ON A HORIZONTAL SURFACE / 
THE SAME IN OUTER SPACE 
HeeeeeDAILY TOTAL RADIATION ON A HORIZONTAL SURFACE 
TOveeeDATLY AVERAGE TEMPERATURE 
XLATs «LATITUDE OF THE LOCATION 
RHOeeeREFLECTIVITY OF THE GROUND AROUND THE WINDOW 
IGeeeeGLAZING INDEX 
1. SINGLE GLAZE 
2. DOUBLE GLAZE 
SCeeeeSHADING COEFFICIENT 
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DIMENSTON LDAY(12)/310289310369310300319319300319°30031/ 

DIMENSTON SMONTH(1209) eSYEAR( 9) »SUMMER(9) 

DIMENSION XDEC (12) 7-19.51 9-10 .280 020011 e560 20 14923027 02N eo 269120 NF 
Kp eSTHM1LN G79 M1965 RIK 23027/ 

DIMENSION R(12)7160391.020701-00579 698759 97270 6 9670 6 96920 9785» 
*—o994591.0133910026791-60327/ 

REAL US(12)/101391013910139101391 006910601 6N601e0691 0601 el3eleI3 
*9 1613/7 °SUM(9) 0H(12) »2KT (12) e TO(12) ep TOT (12) 

DIMENSTON WAZ(B) 0004569900 013556 91806. 9-135 6 9-90 0 P45 oS 

REAL UD (12)72659 66590 665? 6650 0610 e619 0619 e610 0610 0659 0659 0 65/ 


REAL ORTENT(209)/* SOUTH Vet SOUTAWEST,) Men WEST #»* NO 
*RTHWEST %»* NORTH *¢* NORTHEAST ‘%o* EAST %»* SOUTHEAS 
xT * HORIZ. V7, 


REAL TOWN(4) #XSATGS(9) eXSATGD(9) op ZIT(2409) sMLITE(12) #XSG6(9) »xDG(9) 
REAL RST/442.1/9PI/3.1415927/ 9 LAT 
DIMENSION R12) 701426 ol 4s 01567 0189 61967 02059 o 2079 o 201 o 1779 0160 01499. 142/S 
*490,1420/ 
REAL DNI (24) »pAST (2409) oRSI (2409) eSUMD (9) »SUMR(9) 
M=0 
XLAX=AINT(XLAT) 
LAT=(XLAX+ (XLAT-XLAX)/066) *PI/180. 
LAX=INT(XLAT) 
MINUTE=(XLAT=XLAX) *100 
30 FORMAT(*® LOCATION "UAGs 9 LATITUNE'I4*e* NEGREES'I4+e* MINUTES'/) 
CALL DIVOVF (0°00) 
DO 1 N=1012 
RD=AINT(XDEC(N) ) 
DEC=(RDB+(XDEC(N) -RD) /0.6) *PI/180. 
COSWS==TAN (LAT) * TAN (DEC) 
IF (COSWS eGT ele cOReCOSWSeLTerle) RETURN 
IF(NeEQel) WRITE(6041) 
IF(NceEQe1) WRITE(6930) TOWNe LAX e MINUTE 
WS=ACOS (COSWS) 
Y=COS(LAT) *COS (DEC) *SIN( WS) 4WS*STN(LAT) *SIN(DEC) 
HO=24..*R(N) *Y/PI*RST 
TWS=WS*12/PI 
SUNRIZ=12.-ARS(TWS) 
SUNSET=12.+ABS(TWS) 
COSLD=COS (LAT) *COS (DEC) 
SINLD=SIN(LAT) *SIN (DEC) 
S=0. 
DO 500 L=1939 
WW=WS*L/40 6 
CZE=COSLND*COS (WW) +STNLD 
PAR==B(N) /CZE 
APA=ABS (PAR) 
IFCAPA.STe80-) GO TO S501 
ANS=EXP (PAR) *CZE 
GO TO 502 
501 ANS=06 
502 S=ANS+tS 
500 CONTINUF. 
ANSO=E XP (=B(N) 7 £COSLD+SINLD) )* (COSLN+SINLD) /2. 
AI=WS/ 40.*(ANSO+S) 
ZKD=ZD(2KT(N)) 
HO=24./PI*R(N) *RST* (COSLD*SIN (WS) +WS*SINLD) 
THH=HO*ZKT (N) 
DHH=HO* (ZKT(N) =ZKD) 
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23 
24 


25 


26 
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RHH=HO*ZKD 

A=DHH/ (24./PI*AI) 
FAC=A/Z2KT(N) 

DO 2 K=1°9 

DO 2 I=19°24 

DNI(1I)=0. 

ASI(I*#K)=0- 

RSI(IeK)=0. 

ZIT(IeK)=0~6 
DLITE(N)=2-*ABS(TwS) 

DO 3 I=1°24 

TIME=I=-1. 
WT=ABS(12.-TIME) 
W=WT*PI/126 

IF (TIME-SUNRIZ) 30304 

IF (TIME-SUNSET) 50303 
COSZ=SINLD+COSLD*COS (w) 
COSW=COS (DEC) *SIN(W) 
COSS=SQRT (1 .—-COSW*COSW=COSZ*C0SZ) 
V=TAN(DEC) /TAN(LAT) 
TEST=COS (WwW) =V 

IF(TEST) 99908 
COSS==-COSS 
ALT=ASIN(COSZ2) 
AZM=ASIN(COSW/COS (ALT) ) 
IF(COSS) 23924024 
AZM=PI-AZM 

IF (TIME -LTe120¢) AZM=—AZ2M 
AZMP=A2M*180./PI 

DO 10 K=109 

IF(KeNE»9) GO TO 25 
COSTH=C0SZ 

GO TO 26 

CONTINUE 
SAZM=AZM-=WAZ(K)/180.*PI 
SAZMP=SAZM*180-/PI 
ALTP=ALT*180-/PI 
COSTH=COS(SAZM) *COS (ALT) 
CONTINUE 

IF (COSTH-LE-0-) COSTH=0. 
THP=ACOS (COSTH) *180./PI 
PAR2==-B(N) /COSZ 
AP2=ABS(PAR2) 
IF(AP2.G6T-80.) GO TO 3 
DNI(I)=A*EXP(PAR2) 
IF(DNI( I) LE«O-.) DNI(I)=0-6 
DHI=DNI (TI) *COSZ 
IF(DHI.LE-0-) DHI=0. 

ASI (I eK) =DNI(1I) *COSTH 

IF (AST (IeK) LEO.) ASI(IeK)=0. 


RR=P1I/24e* (COS (W)-COS (WS) )/(SIN(WS) =WS*COS (WS) ) 


IF(CRReLTe0-) RR=0. 

RHI=RHH*RR 
RSI(T*#K)=(RHI+(RHI+DHI) *RHO) /26 
ZIT(IeKI=ASI (I eK) +RSI (I eK) 
CONTINUE 
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40 
4e 


31 


506 


14 


505 


32 
35 


CONTINUE 

IF (McEQe1) GO TO 42 

DO 40 T=4019093 

IF(NeEQeI) WRITE(6041) 

IF(NeCEQe I) WRITE(6930) TOWNeLAXe MINUTE 

IF(TeEQe4) WRITE (67101) 

CONTINUE 

IF(M.EQe0) GO TO 43 

IF(NceEQe1) GO TO 43 

WRITE (6041) 

WRITE (6030) TOWNeLAX »MINUTE 

WRITE (60101) 

FORMAT (1H1) 

FORMAT (5F10.3°15e7F10.3) 

WRITE (6031) NeoHOrTO(N) ep TOTIN) » ZKT(N) » 2KD 

FORMAT(/* MONTH* I4o? HO*F7.0e? TON"F7.0e? TOT*F7 0° 
*? KT°FI030 : 
*x? KD°FL1O0¢63//" ORIENTATION’ 31Xe "HOURLY SOLAR RADIATION? 31X0 "DAI 
*LY TOTAL SAG/DAYTIME SAG/ALL DAY? 
* 12X04 5 6 7, 8 9 10 1 
*1 le 13 14 US 16 17 18 19 20 RADIATION SATG 
*S SATGD SATGS SATGD*/11X083(%=")) 

DO 12 K=1909 

SUMN=0. 

SUMD(K)=0-6 

SUMR(K)=9. 

SUM(K)=06 

DO 14 J=1024 

SUMN=SUMN+DNTI (TI) 

SUMD(K) =SUMD (K) +ASTCI¢K) 

SUMR(K) =SUMR(K) +RSI (I eK) 

SUM(K)=SUM(K) +Z1IT(I eK) 

IF (MceEQe Le ANDeKeEQe1) WRITE (60505) (ONI (I) » I=5921) 9 SUMN 
FORMAT (* NORM.DIR. °F4e0016F5-00F10.0) 

SMONTH (Ne K) =SUM(K) *e87*LDAY(N) 412000. 
SATGS=TO(N) +SUM(K) ZOLITE(N) JUS (N) *0287 
SATGD=TO(N) #SUM(K) JDOLITE (N) *SC/UD(N) *0.87 

XSATGS(K)=SATGS 

XSATGD (K)=SATGD 

XSG(K) =TOT(N) +SUM(K) /DLITE(N) /US(N) *0 287 

XDG(K) =TOT(N) #SUM(K) /DLITE(N) *SC/UD(N) *0.87 

IF (M.EQ-0) 
*WRITE (6032) (ORIENT (J0K) 9 J=102) o (ZIT(IeK) 0 L=5021) oSUM(K) » XSATGS(K) » 
*XSATGD(K) »XSG(K) » XDG(K) 

IF (MeEQe1) 
*WRITE (6935) (ORIENT (UeK) »J=102) o (ZIT(C IK) ep I=5 021) »SUM(K) »XSATGS(K) » 
*XSATGD(K) »XSG(K) »XDG(K) 

FORMAT (A6rA4PFH eC OPLEF5 e090 FLOe00 4K 4F6.0) 

FORMAT (/2A6/3X0 "TOTAL "17F5e00F10¢004Xe4FE.0) 

IF(MeEQ-0) GO TO 12 

WRITE (69503) (ASI (IeK) 0 1=50 21) » SUMD(K) 

FORMAT (3X0 *DIRECT°17F5.00F10-0) 

WRITE (60504) (RSI (I¢K) e 1=5921) »SUMR(K) 

FORMAT (3X0 *DIFFUS°17F5-0°F10.0) 

CONTINUE 

CONTINUE 

RETURN 

END 
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FUNCTION ZD(2T) 
DIMENSION ZKT(6)/630 040 oSeco6r ele oT 5/ 


DIMENSION ZKD(6)/.1790 61830 01889 61749 01499 .125/ 


OG £7 


IF(ZT-0-3) 1eie2 
ZD=.179 
GO TO 10 


IF(2Z2T-0-75) 39304 


ZD=.125 

GO TO 10 

DO 20 J=296 
T1I=ZT—-2KT (JU-1) 
T2=ZT=-2KT (J) 
TEST=T1*T2 
IF(TEST) 506920 
Y1=ZKD(J=-1) 
Y2=Z2KD (J) 


ZD=Y1+(V2-Y1) *(ZT=ZKT (J-1) SC Z2KT (J) =2KT (J=1) ) 


GO TO 20 
IF(T1) 89998 
ZD=ZKD(J~1) 
GO TO'20 
ZD=ZKD (J) 
CONTINUE 
RETURN 

END 
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